en ■ 

o. 

Oh. 

<: 



O 
6 



> 

o 
o 

^' 

o 
m 



1 



Mon. Not. R. Astron. Soc. 000,[T]tT7](2013) Printed 29 April 2013 (MN KTeX style file v2.2) 

On the internal pollution mechanisms in the globular cluster NGC 
6121 (M4): heavy-element abundances and AGB models* 

V. D'Orazi^'^t, S.W. Campbell^ M. Lugaro^ J. C. Lattanzio^, M. Pignatari^ 
E. Carretta"^ 

'Department of Physics and Astronomy, Macquarie University, North Ryde, NSW 2109, Australia 

^Monash Centre for Astrophysics, School of Mathematical Sciences, Building 28, Monash University, VIC 3800, Australia 

^Department of Physics, Basel University, Klingelbergstrasse 82, Basel, 4056, Switzerland 

"* INAF Osservatorio Astronomico di Bologna, via Ranzani 1, 1 50100, Bologna, Italy 

Accepted 2013 April 25. Received 2013 April 24; in original form 2013 January 30 



ABSTRACT 

Globular clusters display significant variations in their light-element content, pointing 
to the existence of a second stellar generation formed from the ejecta of an earlier genera- 
tion. The nature of these internal polluters is still a matter of debate: the two most popular 
scenarios indicate intermediate-mass AGB stars and fast rotating massive stars. Abundances 
determination for some key elements can help distinguish between these competitor candi- 
dates. We present in this paper Y abundances for a sample of 103 red giant branch stars in 
NGC 6121. Within measurement errors, we find that the [ Y/Fe] is constant in this cluster con- 
trary to a recent suggestion. For a subsample of six stars we also find [Rb/Fe] to be constant, 
consistent with previous studies showing no variation in other .^-process elements. We also 
present a new set of stellar yields for intermediate-mass AGB stellar models of 5 and 6 solar 
masses, including heavy element s-process abundances. The uncertainties on the mass-loss 
rate, the mixing-length parameter, and the nuclear reaction rates have a major impact on the 
stellar abundances. Within the IM-AGB pollution scenario, the constant abundance of heavy 
elements inside the cluster requires a marginal s-process efficiency in IM-AGB stars. Such a 
constrain could still be satisfied by the present models assuming a stronger mass-loss rate. The 
uncertainties mentioned above are limiting the predictive power of intermediate-mass AGB 
models. For these reasons, at the moment we are not able to clearly rule out their role as main 
polluters of the second population stars in globular clusters. 

Key words: stars; abundances; Galaxy: globular clusters: individual (NGC 6121); stars: AGB 



1 INTRODUCTION 

In the last decades a large number of observational (both photomet- 
ric and spectroscopic) and theoretical studies have been devoted to 
disentangle the complex nature of globul ar clusters (GCs). With 
a few outstan ding exceptions (a) Centauri. llohnson & Pilachowski 



20 id: M22, iMarino et al 



20081 : M54. ICarrettaetal 



20091: NGC 1851. lYong & Grundahl 



2010bh . GCs are homogeneous in their 



Fe-peak, heavy a- (e.g, Ca and Ti) and trans-iron elements pro- 
duced by the slow and the rapid neutron-capture processes (the s 
and the r processes, e.g., Y , Zr, La, Eu; see Gratton et al. 2004, 2012 
for extensive reviews, and Armoskv et al.ll 19941 James et al.ll2004l : 
ISmithll2008l : iD'Orazi et al.ll2OI0l : lRoederedl201 ll for the analysis 
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focused on the heavy elements in GCs). On the other hand, GCs 
exhibit significant star- to-star varia tions in their l ight-element con- 
tent, C, N, O, F, Na, Mg, and Al l[Cohei]|l999l : llvans et alj[T999l : 
ISmith et alj[2005l : lKavser et alj|2008l : Carretta et al. 2009a,b). De- 
pletions in C, O, and Mg abundances appear together with enhance- 
ments in N, Na, and Al (the light-element anticorrelations). These 
chemical featu res, shared by both unevolved (m ain-sequence and 
subgiant, e.g.. iGratton et al.l 120011 : iRamirez & C ohen 20021) and 
evolved st ars (from the red giant branch -RGB- t o the horizon- 
tal branch, iMarino et ai]|201ll : iGratton et aUbOllI) . indicate that 
(at least) two different stellar generations are currently coexisting 
in GCs. The first generation (FG) stars are characterised by high 
O (C and Mg) and low Na (N and Al) abundances and display 
the same chemical composition of field stars at the correspond- 
ing GC metallicity. The second generation (SG) stars present in- 
stead a N/Na/Al-rich (C/0/Mg-poor) pattern. This indicates that 
they were forged from the ejecta of a fraction of FG stars inside 
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the GC JD'Ercole et"Zll2008l: ICarretta et al.ll2010ch . which must 
have experienced H burning at high temperature (T >30 MK). 
These stars, more massive than those we presently observe in GCs 
{M ~0.8 Mq), had time to evolve and internal temperatures high 
enough to activate in their interiors the CNO, NeNa, and MgAl 
cycles to enh ance N, Na, and Al at the expen se of C, O, Ne, 
and Mg (e.g.. lDenisenkov & Denisenkovalll989l ). The astrophysi- 
cal site where this occurred is still under discussion, with two main 
candidates: intermediate-mass asymptotic giant branch stars (IM- 
AGB, 4Me £ M < 8 Mq) undergoing Hot Botto m Burning (HBB, 
IVenturaetaLll200lll2002l:lD'Ercoleetal.||2010h. or fast rotating 
massive sta rs (FRMS, iDecressin et alj |2007|. 120081 : iKrause et all 
120121120130 . Probing the nature of the internal polluters in GCs is 
fundamental to this field of research because of the strong implica- 
tions related to cluster formation and early evolutionary properties. 

The Li abundances provide us a powerful tool to disentangle 
between these two different scenarios due to the fragile nature of 
this element. At the high temperatures where the CNO cycle oc- 
curs, it is expected that all Li is destroyed (Li starts burning above 
T» 2.5 X lO*" K). Interestingly, while FRMS destroy Li, IM -AGB 
stars can also produce it via the ICameron & Fowlen (119711) mech- 
anism. As a consequence, the simultaneous abundance determina- 
tion of Li and of the elements affected by proton captures (e.g., 
O, Na, Mg, Al, hereafter p-capture elements) may supply quite 
stringent observational con s traints to the origin of the polluters 
JD'Antona & VenturallJOOsI : IVentura & D'Antonall2010t) . In other 
words, if there is no Li production within the polluters, we should 
expect a positive correlation between O and Li and a Li-Na anticor - 
relation jpasguini et aLll2005l ; lEind et alj|2009l : IShen et al.ll2010l) . 
Similarly, ascertaining the abundances of F within a GC is critical 
because the prod uction/destruction of F is heavily depen dent on the 
stellar mass (see lSmith et al.ll2005l ; lD'Orazi et al.ll2013l for detailed 
discussions on this topic). Along with Li and p-capture elements, 
surveys of the i-process elements in GCs establish an extremely 
effective tracer, because they can deliver further information on 
the mass range of the polluter. In the Solar System, three differ- 
ent .s-process components have been identified to contribute to the 
abundances above Fe. (0 The weak i-process component occurs 
in massive stars (M >8 Mq), during the conve ctive core He burn- 
ing phase and the subsequent C she ll burning JRaiteri et al.ll 1991c 
iThe et alj2007l : |Pignatari et al.l2O10l) . In the Solar System material 
the weak component accounts for a major fraction of the ^-process 
isotop es between Fe and Sr (60<A<90, see e.g., iKappeler et alj 
uOllf) . («') The main 5-process component takes place in thermally- 
pulsing low- and intermediate-mass AGB stars and it is respon- 
sible, through the ''C(a,n)'*0 and ^^Ne(a,n)^^Mg reactions, for 
the production of the i-process elements between Sr-Y-Zr and Pb 
JBusso et alJI 19991) . Finally, (iii) AGB stars with low initial metal- 
licity contribute to about 50% of the solar -"^Pb and to most of 
the solar 5-process Bi, defined as the strong ^-process component 
JGallinoetaTIl 19981) . 

Neither of the candidate polluters seem able to successfully 
reproduce all the observe d features in GCs (e.g., iFenner et alj 
1 20041 : iKarakas et alj 120061 ). though a growing body of evidence 
seems to conve r ge towards IM-AGB star s as the main polluters 
JRenzinil l2008l : Id' Orazi & Marind l2010t) . On the other hand, 
based on observations of .s-process elements, Villanova & Geisler 
(2011, hereafter VGll) proposed th at in the mildly metal-poor 
GC M4 (rFe/H1=-l. 16. ikarrisl [19961 -2010 update) the polluters 
were massive main-sequence stars. By analysing FLAMES-Giraffe 
(R~ 20,000) spectra for a sample of 23 RGB stars, complemented 
by another 23 stars from lMarino et alj (2008, hereafter Ma08) ob- 



served with UVES (R~45,000), they derived abundances for sev- 
eral key elements, such as Li, Na, the sum C+N+O, Y, Zr, and 
Ba. In addition to confirming previous results on the light element 
inhomogeneities, they found that, within the observational uncer- 
tainties, all the heavy elements do not vary, with the exception of 

Y They found a difference in the [Y/Fe] ratio of ~0.2 dex between 
EG and SG stars and invoked the weak i-process from massive stars 
as responsible for the Y enhancement in the second population. 
However, the weak .s-process is mostly made in massive stars at 
solar-like metallicity. Given the secondary nature, the efficiency of 
the classic .s-process in massive stars is reduced as metallicity de- 
creases, and its contrib ution to the Y inventory is negligible at the 
metallicity of M4 (e.g.. |Prantzos et al.lll990l : lRaiteri et al.l 19931) . 

Furthermore, in the FRMS pollution scenario massive stars 
can contribute to the light-element inter-cluster enrichment only 
when they are on the main sequence. They do this via a slow me- 
chanical wind, delivering to the interstellar medium the ashes of 
the H burning produced via the CNO cycle. On the other hand, 
i-process elements ai^e produced in the later stages of stellar evo- 
lution (during the He core and C shell burning phases), and ex- 
pelled into the interstellar medium during the core-collapse su- 
pernova explosion (SNII). Theory suggests that all but perhaps 
the most massi ve GCs (> lO^Mp, like a i Cen) caimot retain their 
SN ejecta (e.g.. lRecchi & Danzigenl2005l) . In fact, no variations in 
the Q--ele ments are detected in M4, as well as in other "archety- 
pal" GCs JCarretta et alJuOlOah . Indeed, the FRMS scenario pre- 
dicts no variations in th e ^-process element abundances within GCs 
JDecressin et al J 120071) . The only possible contribution via stellar 
winds to i'-process elements from massive stars is from stars that 
become type C Wolf-Rayet stars, where the envelope and the H- 
shell layers have been previously lost. However, it is only in the 
advanced stages of the evolution of these stars that the He she ll 
may show some .5-process enrichment (e.g.. lRauscher et al.ll2002h . 

A further constraint on the possible polluters is given by 
the observations of Li. D'Orazi & Maiino (2010, DM10) deter- 
mined the Li content in M4 for a sample of 104 RGB stars, of 
which 32 were below the RGB bump luminosity. They did not de- 
tect any Li-Na anticoffelation, with EG and SG stars sharing the 
same Li abundance. The average values are A(Li)=1.34±0.04 and 
A(Li)=1.38±0.04, respectively for Na-poor and Na-rich stars. This 
implies that Li should not have been destroyed, ruling out FRMS 
(since they destroy Li) and favouring IM-AGB stars, which can 
synthe sise it. Similar results were obtained by iMucciarelli et alj 
1 1201 ih . who focused on Li abundances from the main-sequence 
up to the RGB and pointed out that there is no variation in the Li 
content across the different stellar generations. The same finding is 
also reported by VGll, who found identical Li values for N-rich 
(A(Li)=0.97+0.04) and N-poor (A(Li)=0.97+0.03) stars. 

Finally, whatever i-process donor we consider, it would be 
difficult to explain a chemical pattern bearing an enhancement in 

Y without a similar signature in the neighbour elements Sr and Zr. 
To shed light on this rather confusing picture we present Y abun- 
dances for a sample of 103 RGB stars; for a small sub-sample of 
six stars we were also able to gat her Rb abundances, whose only 
study available so far comes from lYong et al.l (l2008ah . This paper 
is organised as follows: in Section|2]we describe the sample and the 
abundance analysis procedure; in Section|3]we illustrate our results 
and compare it to previous studies. In Section|4]we provide a sum- 
mary of the current knowledge of the abundance pattern in M4, as 
unveiled from several years of high-resolution spectroscopic stud- 
ies, present new IM-AGB models, and discuss their strength and 
weakness in reproducing the observed features together with the 
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Table 1. Line list, atomic parameters and abundances for the Sun and Arc- 
turus 



Specie 



A 
(A) 



X 
(eV) 



log gf A(X)o A(X)„Bo. 



Yn 4883.684 

Yn 4900.120 

Rb I 7800.268 



1.083 
1.032 
0.000 



0.07 

-0.09 

0.13 



2.24 
2.24 
2.60 



1.32 
1.42 
1.98 



impact of stellar and nuclear uncertainties. In Section|5]we present 
our conclusions. 



2 SAMPLE AND ABUNDANCE ANALYSIS 

Our sample comprises 103 RGB stars whose stellar parameters, 
metallicity and p-capture elements are presented by Ma08, while 
Li abundances are given in DM10. We refer to Ma08 for details on 
target selection, observations, and data reduction procedures. Here 
we just recall that spectra were obtain ed with FLAMES @VLT 
(mounted at UT2. IPasguini et al.|[20o3) . fiber feeding the UVES 
high-resolution spectrograph (nominal resolution R=47,000). The 
RED standard setup at 580nm was employed, providing a spectral 
coverage from 4760-6840 A; the typical S/N ratios of our targets 
range from 100-200 per pixel. For a small sub-sample of six stars 
we could also exploit UVES spectra acquired by E. Carretta and 
collaborators (ESO program 085.D-0205); the 860nm setup results 
in a spectral coverage from 6600 to 10600 A, allowing us to include 
the Rb I line at 7800A. 

We derived Y and Rb abundances t hrough the sp ectral syn- 
thesis technique, using the co de MOOG JSnedenll 19731 . 2011 ver- 
sion) and the Kurucz ( 119931) grids of stellar atmospheres, with 
the overshooting option switched off. Although the only stable 
isotope of Y (^'Y) has an odd mass number, the level splitting 
is essentially negligible, because of the small spin and magnetic 
momentum of the Y nucleus. Thus, we adopted a single-line 
treatment focusing on the features at 4883.68 A and 4900.12 A. 
The line atomic parameters are listed in Table [T] where we also 
report our abundances for the Sun and Arcturus {a Bootes). 
By adopting reffo=5770 K, log go=4.44, A(Fe i)=7.52 and mi- 
croturbulence ^q=0.9 km s~', we gathered A(Y)0=2.24 from 
both lines, which is in excellent agreeme nt with values from 
lOrevesse erahl ( Il996l) , |Asplund et al. ('2009). as well as with the 
meteoritic estimates by iLodders & Palme ( 2009) . Regarding Arc- 
turus we inferred A(Y)=1.37±0.04, using as input stellar parame- 



ters Teff=4286 K, log g=1.67, [F e/H]=-0.52, and ^=1.74 km s 



dRamirez & Allende Prietdl201 1 ). This estimate agrees very well 
with that derived by ISmith et al.l ( 120001) . who focused on lines in 
the yellow-red spect ral window and retrieved A(Y) = 1.40±0.15 (see 
iD'Orazi et al.ll201ll for more details on the i-process elements in 
Arcturus). An example of spectral syntheses of the Y ii lines for 
the sample star #29848 is shown in Figure[T] while in Figure|2]the 
resulting abundances for the 103 stars are plotted as a function of 
the stellar parameters. 



star #29848 
"A(Y)=1,44,1.64, 1.84 




4683 4884 4885 

Wavelength (A) 




The Rb abundance analysis was carried out by synthesising the 
resonance line at 7800.268 A, due to the blending with the 



star #29848 
A{Y)=1.40,1.60, 1.80 
4898.5 4899 4899 5 4900 4900.5 

Wavelength (A) 

Figure 1. Spectral synthesis for the Y n lines at 4883.68 A and 4900.12 
A for star #29848. 



high excitation {x=(>.\16 eV) Si i line lying on the left wing 
of that feature (at /i=7799.996 A). We took accurate wave- 
lengths and relati ve line strengths of th e hyperfine structure 
components from iLambert & Luckl ( Il976l) . adopting the terres- 
trial isotopic mixture of ^^Rb/*^Rb=3 (e.g., iTomkin & LambertI 
|l99 W). We included CN features from B. Plez (private com- 
munication) which significantly improved the continuum fitting; 
note, however, that the impact on the Rb synthese s is negligible. 
Our solar an alysis results in A(Rb)o=2.60, as in lOrevesse et alj 
( Il996l,ll998h . while for Arcturus we derived [ Rb/H]=-0.62 which 
agrees very well with values publ ished by iTomkin & Lamberli 
(1999, [Rb/H]=- 0.58), ISmith et all (2000, [Rb/H]=-0.52), and 
I Yong et al.l ( 120051 [Rb/H]=-0.55). In Figure|3]we show an example 
of the Rb spectral synthesis for star #28356. 
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Figure 2. Y abundances as a function of stellar parameters (Teff , microturbulence f , and log g). 
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Star #28356 
A(Rb)=l. 48. 1.78,2.08 



7796 7798 7800 7802 

Wavelength (A) 
Figure 3. Spectral synthesis for the Rb i line at 7800 A for star #28356. 



3 RESULTS 



Our results are presented in Table |2] where we report the 
stellar parameters, the Na and O abundances from Ma08 and our 
[Y/Fe] ratios along with the coixesponding internal uncertainties 
(computed as detailed in Section [2j. Considering the whole sam- 
ple of 103 stars, we inferred a mean value of [Y/Fe]=0.436±0.004 
(rms=0.038). Our main result is that the scatter is significantly 
smaller than the internal uncertainties, indicating that there is no 
Y variation among our sample stars. If we divide FG and SG 
stars, according to their [Na/Fe] ratios as done by Ma08 and 
VGll (with FG stars defined by having [Na/Fe]<0.23 dex), we 
obtain [Y/Fe] =0.434+0.006 (rms=0.034) and [Y/Fe]=0.437+0.005 
(rms=0.038) for FG and SG stars, respectively. Our finding implies 
that in M4, while we observe a conspicuous variation in p-capture 
elements, the two different stellar generations share the same Y 
abundances; this is shown in Figure [4] where the [Y/Fe] ratios are 
plotted as a function of [Na/Fe] and [0/Fe] from Ma08. The varia- 
tions among different stars of 0.58 and 0.35 dex for Na and O are 
not accompanied by a change in the Y abundances (see discussion 
in Section|4]l. 



Two kinds of internal errors alfect the abundances derived from 
spectral synthesis, those due to the best fit determination and those 
related to the input stellar parameters. In relation to the best fit de- 
termination, we assumed as a conservative estimate the standard 
deviation from the mean, as given from the two different Y ii lines; 
typical values range from 0.02 to 0.10 dex. The sensitivity of the 
[Y/Fe] ratios to the adopted set of the Teif, log g, [A/H], and ^ pa- 
rameters were then assessed in the usual way, by changing one pa- 
rameter at the time and inspecting the coiTesponding variation on 
the resulting abundances. Due to the strengths of the Y ii lines, the 
microturbulence is by far the dominant contribution to this error, 
resulting in 0.07-0.08 dex. Following the error estimates by Ma08 
(i.e., Areff=40K, Alog ^=0.12, A^=0.06 km s"', and A[A/H]=0.05 
dex), we found that total uncertainty in the [Y/Fe] ratios due to 
stellar parameters range from 0.08 to 0.10 dex. Given their inde- 
pendence, we then summed in quadrature the two kind of errors, 
providing the total internal errors as given in the last column of Ta- 
ble|2] Uncertainties in the Rb abundances were evaluated following 
the same approach; because our analysis is based on only one line 
(the Rb I resonance line at 7947. 6A is too weak to be detected in 
our sample stars), the errors due to the best fitting procedure are 
0.10 dex. The total error is given in Table|3] 



The evidence we collected from our measurements on the con- 
stancy of Y within M4 is in disagreement with results obtained by 
VGll. The 23 RGB stars observed with UVES by VGll are also 
included in our sample. In Figure|5]we compare the Y abundances 
from the two studies: there is a discrepancy between the two es- 
timates, which is larger for the Na-poor (i.e., FG) stars. Focusing 
only on their UVES sample, VGll found [Y/Fe]=0.20±0.03 and 
[Y/Fe]=0.34±0.02 for FG and SG stars, respectively. 



There are no obvious explanations for the origin of this 
difference since we used the same spectra, the same Y ii line 
at 4900A, model atmospheres, and stellar parameters (all from 
Ma08), as well as the same abundance code, employing the spec- 
tral synthesis technique. In Figure |6] we show the direct compar- 
ison of the spectra for two stars with very similar atmospheric 
parameters, and for which VGll have derived a difference in the 
[Y/Fe] ratios of A[Y/Fe]=0.27, being [Y/Fe]=0.17 for #29848 and 
[Y/Fe] =0.44 for #26794. As can be seen from the figure, the Y ii 
lines are of the same strength in both stars and we indeed ob- 
tained [Y/Fe] =0.43±0. II and [Y/Fe]=0.49±0. 10, respectively. Fur- 
thermore, the low [Y/Fe] ratios derived by VGI 1 for FG stars (rang- 
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Table 2. Stellar parameters, [0/Fe], [Na/Fe] ratios (from Ma08) and Y abundances for our 103 sample stars with the related intemal errors. 



Star 


(K) 


logs 


(kms-') 


[Fe/H] 


[0/Fe] 


[Na/Fe] 


[Y/Fe] 


Internal 
enor 


19925 


4050 


1.20 


1.67 


-1.02 


0.28 


0.51 


0.42 


0.11 


28103 


3860 


0.50 


1.62 


-1.08 


0.50 


0.17 


0.44 


0.10 


33414 


4840 


2.51 


1.28 


-1.05 


0.52 


0.21 


0.35 


0.13 


35508 


4780 


2.48 


1.18 


-1.05 


0.25 


0.38 


0.45 


0.10 


29272 


4780 


2.50 


1.26 


-1.11 


0.53 


0.05 


0.37 


0.10 


28797 


4640 


2.35 


1.36 


-1.12 


0.35 


0.44 


0.43 


0.09 


29848 


4780 


2.52 


1.24 


-1.05 


0.54 


0.09 


0.43 


0.11 


5359 


4800 


2.44 


1.28 


-1.03 


0.42 


0.13 


0.40 


0.09 


20766 


4400 


1.80 


1.45 


-1.05 


0.31 


0.53 


0.40 


0.09 


21191 


4270 


1.60 


1.60 


-1.06 


0.34 


0.51 


0.40 


0.10 


21728 


4525 


2.00 


1.42 


-1.06 


0.31 


0.37 


0.40 


0.09 


22089 


4700 


2.28 


1.36 


-1.06 


0.35 


0.50 


0.45 


0.10 


24590 


4850 


2.66 


1.35 


-1.07 


0.50 


0.30 


0.40 


0.09 


25709 


4680 


2.20 


1.38 


-1.13 


0.44 


0.34 


0.49 


0.11 


26471 


4800 


2.40 


1.28 


-1.10 


0.41 


0.09 


0.48 


0.09 


26794 


4800 


2.45 


1.44 


-1.17 




0.36 


0.49 


0.10 


27448 


4310 


1.57 


1.58 


-1.12 


0.51 


0.11 


0.42 


0.11 


28356 


4600 


2.22 


1.53 


-1.14 


0.44 


0.37 


0.50 


0.10 


28707 


4880 


2.74 


1.31 


-1.03 




0.22 


0.41 


0.13 


28847 


4780 


2.40 


1.27 


-1.16 


0.52 


0.08 


0.51 


0.11 


28977 


4680 


2.33 


1.40 


-1.14 


0.39 


0.40 


0.49 


0.09 


29027 


4720 


2.40 


1.41 


-1.10 


0.51 


0.02 


0.40 


0.10 


29065 


4650 


2.10 


1.41 


-1.12 


0.45 


0.17 


0.40 


0.09 


29171 


4880 


2.64 


1.26 


-0.99 


0.22 


0.41 


0.40 


0.10 


29222 


4720 


2.50 


1.36 


-1.04 


0.36 


0.24 


0.45 


0.12 


29282 


4650 


2.30 


1.42 


-1.06 


0.28 


0.42 


0.45 


0.09 


29397 


4600 


1.50 


1.78 


-1.12 


0.24 


0.46 


0.41 


0.14 


29545 


4880 


2.61 


1.18 


-1.06 


0.47 


-0.02 


0.40 


0.12 


29598 


4840 


2.50 


1.40 


-1.06 


0.39 


0.40 


0.51 


0.11 


29693 


4360 


1.10 


1.88 


-1.19 


0.26 


0.50 


0.41 


0.14 


30209 


4880 


2.62 


1.32 


-0.99 


0.42 


-0.05 


0.44 


0.10 


30345 


4850 


2.73 


1.31 


-1.06 


0.32 


0.43 


0.41 


0.11 


30450 


4760 


2.53 


1.35 


-1.00 


0.40 


0.40 


0.46 


0.11 


30452 


4830 


2.56 


1.25 


-1.00 


0.40 


0.06 


0.40 


0.13 


30549 


4830 


2.52 


1.28 


-1.09 


0.45 


0.13 


0.42 


0.11 


30598 


4360 


1.75 


1.47 


-1.07 


0.43 


0.05 


0.44 


0.09 


30653 


4660 


2.30 


1.25 


-1.06 


0.47 


0.15 


0.46 


0.09 


30675 


4830 


2.58 


1.35 


-1.07 


0.36 


0.41 


0.50 


0.10 


30711 


4560 


2.25 


1.46 


-1.01 


0.37 


0.32 


0.41 


0.09 


30719 


4810 


2.65 


1.24 


-1.19 


0.44 


0.42 


0.50 


0.09 


30751 


4430 


1.78 


1.47 


-1.09 


0.40 


0.32 


0.43 


0.09 


30924 


4810 


2.60 


1.28 


-1.09 


0.48 


0.20 


0.45 


0.09 


30933 


4800 


2.63 


1.30 


-1.07 


0.44 


0.44 


0.50 


0.17 


31015 


4800 


2.47 


1.37 


-1.07 


0.50 


0.00 


0.49 


0.11 


31306 


4900 


2.87 


1.33 


-1.11 


0.30 


0.40 


0.48 


0.13 


31376 


4800 


2.59 


1.36 


-1.00 


0.33 


0.43 


0.41 


0.10 


31532 


4770 


2.60 


1.21 


-1.03 


0.47 


0.38 


0.41 


0.11 


31665 


4650 


2.17 


1.34 


-1.07 


0.48 


0.02 


0.42 


0.10 


31803 


4850 


2.60 


1.34 


-1.13 


0.38 


0.25 


0.37 


0.16 


31845 


4700 


2.42 


1.31 


-1.06 


0.27 


0.38 


0.39 


0.09 


32055 


4300 


1.57 


1.52 


-1.12 


0.31 


0.40 


0.38 


0.10 


32121 


4840 


2.58 


1.33 


-0.93 




0.35 


0.46 


0.11 


32151 


4770 


2.43 


1.38 


-1.07 




0.39 


0.44 


0.12 


32317 


4510 


1.88 


1.43 


-1.07 


0.37 


0.37 


0.48 


0.10 


32347 


4640 


2.22 


1.38 


-1.09 


0.40 


0.23 


0.39 


0.09 


32583 


4850 


2.54 


1.34 


-1.08 


0.46 


0.07 


0.38 


0.12 


32627 


4750 


2.42 


1.30 


-1.10 


0.50 


-0.03 


0.46 


0.11 


32700 


4560 


2.12 


1.36 


-1.02 


0.49 


0.06 


0.43 


0.09 


32724 


4850 


2.73 


1.30 


-1.03 


0.51 


0.11 


0.41 


0.10 


32782 


4880 


2.60 


1.24 


-1.05 


0.29 


0.46 


0.44 


0.10 



V. D 'Orazi et al. 



Table 2. Continued 



Star 


Teff 

(K) 


logs 




[Fe/H] 


[O/Fe] 


[Na/Fe] 


[Y/Fe] 


Internal 
error 


32871 


4770 


2.48 


1.24 


-1.01 


0.40 


0.06 


0.47 


0.17 


32874 


4600 


2.04 


1.38 


-1.14 


0.41 


0.11 


0.45 


0.09 


32933 


4430 


1.42 


1.78 


-1.13 


0.48 


0.02 


0.43 


0.13 


32968 


4630 


2.17 


1.30 


-1.13 


0.41 


0.24 


0.45 


0.09 


32988 


4850 


2.63 


1.22 


-1.09 


0.51 


0.02 


0.50 


0.10 


33069 


4940 


3.05 


1.36 


-0.92 


0.39 


0.23 


0.48 


0.09 


33195 


4620 


2.38 


1.26 


-1.03 


0.48 


0.09 


0.50 


0.11 


33617 


4720 


2.35 


1.26 


-1.09 


0.29 


0.37 


0.50 


0.09 


33629 


4930 


2.80 


1.33 


-0.98 


0.40 


0.31 


0.39 


0.09 


33683 


4800 


2.57 


1.18 


-1.05 


0.42 


0.00 


0.48 


0.09 


33788 


4700 


2.37 


1.33 


-1.02 


0.30 


0.37 


0.44 


0.09 


33900 


4770 


2.48 


1.27 


-1.06 


0.30 


0.47 


0.45 


0.09 


33946 


4800 


2.62 


1.15 


-1.03 


0.23 


0.32 


0.40 


0.10 


34006 


4320 


1.67 


1.61 


-1.06 


0.25 


0.44 


0.43 


0.11 


34130 


4550 


2.08 


1.40 


-1.09 


0.33 


0.43 


0.40 


0.09 


34167 


4950 


2.60 


1.40 


-1.10 




0.13 


0.42 


0.13 


34240 


4470 


1.95 


1.41 


-1.10 


0.47 


0.08 


0.42 


0.10 


34502 


4860 


2.70 


1.33 


-1.08 


0.35 


0.34 


0.48 


0.09 


34726 


4600 


2.24 


1.35 


-1.01 


0.37 


0.42 


0.43 


0.09 


35022 


4850 


2.51 


1.36 


-1.08 


0.30 


0.41 


0.39 


0.10 


35061 


4860 


2.67 


1.23 


-0.99 


0.50 


0.06 


0.48 


0.10 


35455 


4600 


2.10 


1.29 


-1.06 


0.45 


0.01 


0.47 


0.09 


35487 


4850 


2.67 


1.24 


-1.00 


0.36 


0.30 


0.47 


0.10 


35571 


4880 


2.79 


1.10 


-0.99 


0.45 


-0.02 


0.45 


0.10 


35627 


4830 


2.37 


1.20 


-1.08 




0.04 


0.42 


0.09 


35688 


4720 


2.25 


1.33 


-1.11 


0.21 


0.40 


0.42 


0.09 


35774 


4450 


1.92 


1.44 


-1.10 


0.43 


0.24 


0.42 


0.09 


36215 


4300 


1.59 


1.53 


-1.11 


0.48 


0.18 


0.43 


0.09 


36356 


4820 


2.66 


1.26 


-1.05 


0.30 


0.26 


0.49 


0.09 


36929 


4820 


2.55 


1.28 


-1.03 


0.34 


0.45 


0.44 


0.10 


36942 


4800 


2.66 


1.23 


-0.98 


0.55 


0.04 


0.43 


0.13 


37215 


4790 


2.50 


1.21 


-1.11 


0.45 


0.25 


0.41 


0.13 


38075 


4800 


2.54 


1.25 


-1.07 




0.42 


0.49 


0.13 


38383 


4590 


1.87 


1.45 


-1.10 


0.39 


0.11 


0.37 


0.13 


38896 


4760 


2.53 


1.31 


-1.02 


0.31 


0.43 


0.43 


0.09 


42490 


4570 


2.08 


1.41 


-1.07 


0.29 


0.38 


0.43 


0.09 


42620 


4600 


2.05 


1.37 


-1.09 


0.22 


0.48 


0.39 


0.11 


43370 


4920 


2.80 


1.32 


-1.05 


0.20 


0.34 


0.41 


0.11 


44243 


4860 


2.80 


1.17 


-1.04 




0.41 


0.47 


0.10 


44595 


4750 


2.40 


1.40 


-1.07 


0.20 


0.44 


0.40 


0.09 


44616 


4620 


2.20 


1.44 


-1.04 


0.28 


0.44 


0.40 


0.10 


45163 


4770 


2.40 


1.26 


-1.10 




0.46 


0.45 


0.12 


45895 


4720 


2.25 


1.34 


-1.04 


0.43 


0.10 


0.43 


0.09 



ing from 0.07 to 0.29) does not fit the global picture of the neutron- 
capture element abundances in M4. In fact, it is well assessed from 
several studies that M4 is characterised by an intrinsically high .s- 
process element content, compared t o other GCs (like e.g., i ts twin 
M5, see lBrown & Wallersteinlll 9921; ll vans et al.lll999ll200ll : Yong 



et al. 2008a.b: Ma08: ISmithll2008l : lD'Orazi et aU201(]h 

On t he other hand, five o f the stars in our sample are in com- 
mon w ith lYong et al.l j2008ah . The average [Y/Fe] in lYong et alj 
( l2008ah and in our study is 0.69±0.02 and 0.42±0.01 dex, respec- 
tively, with a difference of 0.27 dex. Such a difference is mainly 
due to a [Fe/H] offset of 0.18±0.02 dex, whereas the average [Y/H] 
values differ only by 0.08±0.02 dex, within the observational un- 
certainties in the adopted stellar parameters. For instance, T^g val- 



ues of Yong and collaborators are on average cooler than Ma08 and 
log g are systematically lower by -0.43±0.07. 



For 6 stars we derived an average value of [Rb/Fe]=0.34±0.01. 
Although based on a quite limited sample, the Rb a bundance is 
constan t within M4, confirming previous findings by lYong et alj 
( 120080) . One of our star is in common with that study: 
#36215 (L3624 in Yong et al.'s sample) for which we derived 
[Rb/H]=-0.79 to be compared with [Rb/H]=-0.89 by Yong et al. 
In general there is a good agreement between the two estimates: our 
mean value is [Rb/H]=-0.77±0.04 (rms=0.09), to be compared to 
[Rb/H]=-0.84±0.03 (rms=0.09) (see FigurelT}. 
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Figure 4. [Y/Fe] (this study) vs [Na/Fe] and [O/Fe] (from Ma08) for our sample stars. 
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Figure 6. Comparison of the two Y n lines for two stars of almost identical stellar parameters (#29848, solid line; #26794, dotted line). 



Table 3. Rb abundances for a sub-sample of six stars (see text) 



Star 


[Rb/Fe] 


Internal 
error 


28356 


0.35 


0.10 


29693 


0.27 


0.12 


30711 


0.35 


0.10 


30751 


0.36 


0.11 


32317 


0.34 


0.10 


36215 


0.32 


0.12 



4 DISCUSSION 

M4 is perhaps the most extensive ly studied GC . Mainly thanks to 
its close distance (Rsu„=2.2 kpc, lHarrislll99a) . a wealth of spec- 
troscopic studies have been accomplished, aimed at deriving in- 



formation on its abundance pattern for species from Li up to Pb. 
This cluster represents an excellent example of what we deem as an 
archetypal GC, because it presents substantial changes only in p- 
capture elements, while heavy elements do not show any variations 
beyond what is expected from observational errors. A further con- 
firmation of the "standard" nature of M4 can also be found in the 
recent survey by Carretta et al. (2012a,b), who performed a careful 
analysis of the Al content in three GCs. In NGC 6752 and 47 Tuc 
(NGC 104) at least three episodes of star formation are required to 
account for the chemical pattern. In M4 the classical picture of a 
FG plus a SG formed from the ejecta (diluted to a certain extent) of 
FG stars is able to reproduce the observed trends. It is worth men- 
tioning, in this context, that the cluster is however characterised 
by an intrinsically high level of .?-process elements, when com- 
pared to other clusters like M5 (which has a very sim ilar metallicity, 
[Fe/H]« -1.2 dex; see e.g.. llvans et aljl999l . [200ll) . Several studies 
have been devoted to assess this very peculiar heavy element pat- 
tern, indicating that the proto-cloud from which M4 formed must 
have b een enriched by a n higher concentration of .?-process ele- 
ments JYong et alj2008ah . 

In the following Section 143] we briefly overview the current 
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Figure 5. Compaiison between our Y abundances and those by VGll for 
FG (empty squares) and SO stars (filled squares). 
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Figure 7. [Rb/H] ratios as a function of temperature for our stars (squares) 
and for the 12 giants from Yong et al. (2008a, starred symbols) 



knowledge of the chemical composition of this cluster. In Sec- 
tion 14.21 we present a new set of intermediate-mass AGB models 
and discuss their comparison with the observational data to chal- 
lenge the IM-AGB pollution scenario. 



4.1 The chemical abundance pattern of M4 

Light elements 

The most recent abundance determinations for the key elements Na, 



O, Mg, Al were presented by Ma08, Carretta et al. (2009a,b) and 
VGl 1. Ma08 analysed a sample of 104 giants from high-resolution 
UVES spectra and detected a well defined Na-O anticorrelation 
(left panel of Figure [8}, with [0/Fe] ranging from 0.20 to 0.55 
(A=0.35 dex) and [Na/Fe] ranging from -0.05 to 0.53 dex (A=0.58 
dex). Very similar results were obtained by VGll, who confirmed 
0.10<[O/Fe]<0.53 and -0.11<[Na/Fe]<0.57, while slightly larger 
variations in the Na-0 plane were measured by ICarretta et al.l 
( 1200%) . Based on the intermediate-resolution FLAMES-GirafFe 
(R~20000) spectra of 88 giants, and a diff^erent temperature scale 
based on photometry, they found the [O/Fe] ratio varying from 
-0.15 to 0.45 (A=0.60 dex) and [Na/Fe] from -0.16 to 0.70 dex 
(i.e., A=0.86 dex). While all the three works detected the presence 
of a very clear Na-O anticorrelation, there is no consensus on Al. 
Ma08 found a variation in [Al/Fe] (at ~0.4 dex), positively c orre- 
lated with Na, as also previously detected bv llvans et al.l ( ll999f) . but 
no significant (anti)correl ation with the Mg ab undances. At vari- 
ance with that study, both lCarretta et al.l (l2009ah and VGl 1 did not 
detect any significant variation in the Al content of the cluster (right 
panel of Figure[8}. Carretta et al. presented Mg and Al abundances 
for a sample of 14 RGB stars from UVES spectra and obtained con- 
stant values for [Mg/Fe]=0.55±0.01 and [Al/Fe]=0.60±0.01; anal- 
ogously VGl 1 found no hint of Al variation, with FG and SG stars 
showing [Al/Fe] =0.5 1+0.04 and [Al/Fe]=0.53±0.02, respectively. 
On the other hand, Carretta et al. (2012b) presented Al abundances 
for 83 RGB stars derived from the strong doublet at 8772-8773 A, 
properly taking into account the CN features populating that spec- 
tral region. They detected a relatively small variation in Al, posi- 
tively correlated with Na, and confirmed the constancy of the Mg 
abundances ([Mg/Fe]=+0.541±0.005 dex). In general, beyond the 
debated question of the existence of an Al variation, all these stud- 
ies agree that the cluster does not show the Mg-Al anticorrelation. 

As previously mentioned in Section [T] Li a bundances in M4 
were derived from three different studies: DMlO. lMucciarelli et alj 
(2011), and VGll. Despite the presence of systematic offsets be- 
tween the different studies (Li is critically dependent on the adopted 
effective temperature scale), the main result is that there is no Li- 
Na anticorrelation and no Li-O correlation (which would be ex- 
pected in the case of Li destruction within the polluters), since 
both FG and SG stars show the same Li content. While DM10 
and VGll analysed the Li content in RGB stars (bo t h belo w and 
above the RGB bump luminosity), iMucciarelli et al.l ( 120111 ) gath- 
ered Li abundances for a sample of dwarf members, and report 
A(Li)=2.30±0.02 dex (rms=0. 10), consistent with the Spite plateau 
defined by halo dwarfs. The fact that FG and SG stars exhibit the 
same Li calls for yields from the polluters to be at roughly the Spite 
plateau level and led DM10 to speculate that relatively low-mass 
polluters (a4 Mq) were at work in M4. Further support to this sug- 
gestion comes from the modest Na-O anticoiTelation (with the ab- 
sence of extre mely 0-poor stars, as compared to other GCs, such 
as NGC 2808, lBragaglia et alJlIoToh and from the lack of the Mg- 
Al anti-correlations (see, however, the next Section for comparison 
between models and observations). 

Concerning the C and N abundances, the cluster is known to 
exhibit a very well defined bi modal distribut i on in the CN band 
strengths, as first discovered by iNorrisI (119811) . Ilvans et al.l d 1 9991) 
found that variations in CN are not random but they are positively 
con'elated with Na and Al (with the CN-strong stars showing also 
higher abundances of Na and Al) and anticorrelated with O. A 
similar result was later presented by Ma08 who, adopting the CN 
strength index S(3839) as given by Smith & Briley (2005), showed 
the (anti)correlation of CN index with (0)Na (see Figure 10 of that 
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Figure 8. Na-O anticorrelation (left panel) and [Al/Fe] as a function of [Mg/Fe] (right panel). Data are from Ma08 (filled circles); Carretta et al. (2009a,b, 
empty squares); VGll (empty triangles). The relatively small oft'sets between the dift'erent studies can provide a conservative estimate of the measurement 
uncertainties. 



paper). The sum of C+N+0 has been f ound to be constant, within 
the uncertainties, by llvans et al.l ( 119991) and VGll, who obtained 
A(C+N+O)=8.24±0.03 and A(C+N+O)=8.16±0.02, respectively. 
We note, in passing, that the constancy of the C+N+O sum in GCs 
is currently debated. A most striking case is NGC 1851, for which 
Yong (2011) detected an huge spread of about 1.6 de>iJ, while 
Ivillanova et al.l ( 120 lOl) claimed a constant value. 

Finall y, the He conte n t has been measured for six blue HB 
stars from Ivillanova et al.l ( 120 12h . who found these stars all en- 
hanced in He by 0.04 dex with respect to the primordial He abun- 
dance (i.e., Y=0.24-0.25). The same stars are also Na-rich (O- 
poor), as expected acc ording to the multip le population scenario 
(we refer the reader to iGratton et allbOlOl and iMarino et all 1201 ll 
for a discussion on the relationship between light-element varia- 
tions, HB morphology, and the second-parameter problem). 



Heavy elements 



Yong et al. (2008a,b) published a comprehensive abun- 
dance study for a sample of twelve RGB stars in M4, deriv- 
ing abundances for 27 elements, from Si to Hf. These authors 
concluded that there is no change in any of the n-capture ele- 
ments under scrutiny; regarding Y they found a mean value of 
[Y/Fe]=0.69±0.02 (rms=0.05). The FG stars (i.e., [Na/Fe]<0.30) 
exhibit [Y/Fe]=0.67±0.03 (rms=0.05) while Na-rich, SG stars have 
[Y/Fe]=0.70±0.02(ims=0.06). The Pearson's correlation coeffi- 
cient between [Y/Fe] and [Na/Fe] is very small, showing that there 
is no statistically significant relationship between the two. Simi- 
larly, none of the other heavy elements show any significant varia- 
tion, with [Sr/Fe]=0.73±0.04, [Zr/Fe]=0.48±0.03 (from Zr n lines), 
[La/Fe]=0.48±0.03, [Eu/Fe]=0.40±0.03, and [Pb/Fe]=0.30±0.02 

dex. 

Barium a bundances w ere published by llvans et alj 

( Il999l) , Ma08, iD'Orazi et all ( l201(il) . and VGll who found 
[Ba/Fe]=0.60+0.02, [Ba/Fe]=0.41±0.01, [Ba/Fe] =0.50+0.01, and 



^ 0.5 

0) 

\ 

CO 

PQ 




o 



J I I L. 



± 



J I I L. 



J_ 



J I I L. 



-0.5 



0.5 

[Na/Fe] 



Figure 9. [Ba/Fe] vs [Na/Fe]. Symbols are as for Figure |8]for Ma08 and 
VGll, while empty squares and pentagons are for D'Orazi et al. (2010) 
and Ivans et al. (1999), respectively. 



[Ba/Fe]=0.32±0.02, respectively. We show in Figure|9]the [Ba/Fe] 
ratio as function of the [Na/Fe] ratio, reporting measurements 
from these previous studies. Although the Ba abundances are 
characterised by relatively large uncertainties (due to the saturated 
behaviour of the Ba n lines), there is no obvious relationship with 
the Na abundances (nor with other p-capture elements). 

A key element in this context is Rb (Z=37), w hose abundance 
in the solar system material is ~50% s-process dSimmerer et al.l 
l2004h . Theory predicts that AGB stars with M >5Mq overpro- 
duce Rb with respect to solar and the other nearby 5-process el- 
ements (e.g., Y Zr). This is the result of neutron densities reach- 
ing up to 10''* n/cm' and favouring the operation of t he branch- 
ing points on the .y-proce ss path at ^^Kx and ^*Rb dAbia et alj 



mg points on the .y-proce ss path at Kr and Kb dAbia et alj 
' \http : //www.ucolick.org/kraftfest/talks/yong.kraftfest.pdf\ uOOltlvan Raai et al.ll2012h . As previously mentioned in SectionjS] 
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lYong et alj l l2008ah obtained that the Rb content does not show any 
cosmic star-to-star scatter (i.e., [Rb/Fe]=0. 39+0.02, rms=0.07). 
For FG and SG stars, the average ratios are [Rb/Fe]=0.37±0.02 
(rms=0.04) and [Rb/Fe]=0.40±0.02 (rms=0.07), respectively. We 
confirm this previous finding and find a mean abundance of 
[Rb/Fe]=0.34±0.01. 

In summary, all the previous surveys on heavy-element ele- 
ment abundances in M4, with the exception of Y by VGll, agree 
that these elements do not show any intrinsic inhomogeneity, within 
the observational uncertainties. This demands that the polluters 
must account for changes in p-capture elements without affecting 
the n-capture elements (see Section |4~2) . 

Finally, since it is not the goal of this paper to discuss the 
intrinsic enrichment of heavy elements in M4 and to address the 
origin of the first generations of stars responsible for such a pecu- 
liar signature, we refer to previous papers available in the liter ature 
on this topic (e.g., Yong et al. 2008a,b; iKarakas et"aLll2010l') . We 
simply take these abundances of light and heavy elements as ini- 
tial composition and conditions to study the chemical enrichment 
inside M4 (see section 4.2). 



4.2 Intermediate-mass AGB models 

4.2.1 The stellar models 

We have computed a grid of stellar models to investigate the sce- 
nario where the material that went into forming the second gen- 
eration stars is assumed to come from the ejecta of a population 
of IM-AGB stars. To compute the stellar structure models we used 
the Monash version of the Monash-Mount Stromlo evolutionary 
code (MONSTAR see eg. IWood & Zarr"olll98ll ; iLattanziol [19861 : 
iFrost & Lattanzid 1996h includ ing recent updates as described 
bv lCampbell & Lattanzid 12003). Low-temper ature opacities have 
been further updated to those calculated by iLederer & Aringeg 
( 120091) ■ which are variable in C and N as needed to follow the pollu- 
tion of the s tellar surface v ia the third dredge-up (TDU) episodes in 
AGB stars ( lMarigdl2002h . Instantaneous mixing was used in con- 
vective zones and the convective b oundaries were defined u sing 
a search for "convective neutrality" ( IFrost & Lattanziolll99q) . No 
further overshoot was applied beyond this boundary. Initial com- 
position for the models was based on the observed FG abundances 
reported by VGll (see their Table 2). This composition is alpha- 
enhanced and has sub-solar C/Fe. Since the initial [0/Fe] = +0.4, 
the metallicity Z is roughly doubled from that given by [Fe/H] only 
and the present models have initial Z = 0.002. In our stand ard mod- 
els ma ss-loss was included using the empirical formula of lReimersI 
( I1975I) du ring the RGB phase (with 77 = 0.4) and the empirical 
formula of IVassiliadis & Wood ( 119930 during the AGB phase. The 
stellar evolution was followed from the zero-age main sequence 
to the end of the thermally-pulsing AGB phase. In some cases the 
evolution could not be followed past the insta bility caused by the 
peak in Fe opacity in the AGB envelope (see lLau et al.ll2012l for 
a discussion) and final envelope masses between 0.1 M© and 1.2 
Mq remained at the end of the calculations. In these cases when 
calculating the yields we assume that the surface abundances in 
this remaining material stay the same as those in the last computed 
model. 

Models of AGB stars are well known to contain many sig- 
nificant uncertainties. These include, but are not limited to: (1) 
the convection theory, in these models we used the mixing length 
theory (MLT), (2) the mass-loss rates, and (3) the nuclear reac- 
tion rates. Our grid of models explores the effects of these three 



uncertainties. We calculated two 'standard' (STD) models of ini- 
tial masses 5 M© and 6 M© using the input physics described 
above and setting the mixing length parameter a^tr = ImltIHp 
(where Imlt is the mixing length and Hp is the pressure scale- 
height) to 1.75, as calibrated using a model of the Sun. As a first 
test, we increased Umlt to 2.20. Variations of Qmlt to this ex- 
tent are warrante d since Umlt probably changes for diflFerent evo- 
lutionary ph ases (Lvd on et al. 1993; Lebzelter & Wood 20 071) and 
metal l icities dChiefln et al. 1995: Palmieri et al. 2002 : Ferraro et al.l 



120061 : iTrampedach & Stein 20 III) and it is expected to change 
for AGB~stars dSackma nn & Boothrovdl ll99ll: iLebzelter & WoodI 
l2007h . lVentura & D' An tona (2005a) found that the Omlt value cor- 
respo nding to employing the full spectrum of turbulence formalism 
(FST. ICanuto & Mazziteiiilll99lh in massive AGB stars is ~ 2 (al- 
though this does not re produce all the features of FST models). 
iMcSavenev et al.l j2007h used values of autT up to 2.6 to match 
observations of the effective temperature of the envelopes of mas- 
sive AGB stars in the Ma gellanic Clouds. A s a second test, we used 
the mass-loss formula of iBloecken ( Il995h . with r] = 0.02 during 
the AGB phase. This m ass-loss prescription gives h igher mass-loss 
rates than the foiTnula of lVassiliadis & Wood ( Il993h . eifectively re- 
ducing the lifetime of the star. The star consequently suffers fewer 
thermal pulses a nd therefore fewer TDU epis odes, which follow 
each TP (see also lVentura & D'Anton3l2005bl) As a third test, we 
changed both the mass-loss rate to iBloecken ( 119950 and smlt to 
2.2. 

We also performed detailed nucleosynthesis calculations 
using a post-processing c ode, the Monash Stellar Nucleosyn- 
thesis code (MON SOON; ICannod 1 1 9931 : iLattanzio et al] 1 1 9961 : 
iLugaro et al.l l2004r) . that reads as input the background thermal 
structure from the stellar evolution calculations. The code solves 
a network of 2,336 nuclear reactions involving 320 nuclear species 
from n and p up to Bi. The bulk of the reaction rates are from the 
JINA reaclib database as of May 2009. We tested that there are no 
major changes when using the database as of May 2012. During 
the post-processing we made one final test to investigate the ef- 
fects of changing the p-capture rates that affect the production of 
^^Na, as detailed in the last line of Table |4] (model MBR, see also 
IVentura & D'Antonal2005a : 2006). We summarise the 9 computed 
models and their parameters in Table |4] 



4.2.2 Model results 

In the last two columns of Table |4] we list the total number of ther- 
mal pulses and the final core mass for each model. The models 
show little difference in final core mass because the core grows 
very slowly due to the presence of deep TDU episodes, which re- 
moves the core growth due to H burning (i.e., Xjoy ~ I). The mass 
of the core is not affected by higher mass-loss rates since it essen- 
tially reaches its final value after ~ 10 TPs. On the other hand, the 
different models present large differences in the number of TPs. In- 
creasing Omlt caused a minor reduction in number of TPs in the 5 
Mq model while in the 6 M© model the number of TPs was reduced 
by ~ 25%. This is a result of hott er HBB and higher luminosi - 
ties driving an enhanced mass loss dVentura & D'Antondl2005ah . 
In the models w ith increased mass-loss rates, due to the use of the 
iBloeckeg ( 119951) formula, obviously the effect is stronger. In both 
the 5 Mq and 6 Mg models the number of TPs as compared to the 
standard models is significantly lower, reducing from 79 — » 39 TPs 
and 76 ^ 47 TPs, respectively. As expected, the combination of 
a higher a^tT and a higher mass-loss rate leads to an even further 
reduction of the number of TPs (Table 0. 
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Table 4. List of our grid of models. Final Ma,ie is the mass of the remnant white dwarf. 



LABEL 



Initial Mass AGB Mass-loss 



omlt 



Ne Rates 



No. TPs Final M„ 



STD5 
MLT5 
BLK5 
M+B5 



5.0 
5.0 
5.0 

5.0 



Vassiliadis & Wood (1993) 
Vassiliadis & Wood (1993) 
Bloecker (1995) 
moecker (1995) 



1.75 
2.20 
1.75 
2.20 



Standard 



79 
73 
39 
29 



0.91 
0.91 
0.91 
0.91 



STD6 
MLT6 
BLK6 
M+B6 
MBR 



6.0 VassiUadis & Wood (1993) 1.75 

6.0 VassiUadis & Wood (1993) 2.20 

6.0 Bloecker (1995) 1.75 

6.0 Bloecker (1995) 2.20 

6.0 Bloecker (1995) 2.20 



-2Ne(p,y)"Nahigh'' 
& 23Na(p,ff)2'*Mg low" 



"As given by Iliadis et al. (2010) . 



76 

55 
47 
32 
32 



1.00 
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0.99 
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Figure 10. The evolution of the surface C, N, O, Na, Mg, Al, and the sum of 
the C+N+O abundances in the standard 6 M© model (STD6, upper panel) 
and in the model computed using the Bloecker mass-loss rate formula and 
oiMLT = 2.2 (model M+B6, lower panel). Time is offset so / = coincides 
with the start of the AGB for both models. We note that the Al in the plot 
is mostly ^*A1 which, due to its relatively short half-life, was assumed to 
decay completely to ^^Mg and was thus added to the yields of Mg presented 
in Figs.fTilandfTsl 



In Figures [Tolandfm we compare the evolution of the surface 
abundances along the AGB in the standard 6 Mq model (STD6) 
and in the 6 M© model computed using Bl oecker mass-loss and 
ttMLT = 2.2 (M-I-B6). Using the formula of lBloecked ( 1 19951) and 
Omlt = 2.2 severely truncates the AGB phase, resulting in fewer 
TDU episodes and less material carried into the envelope. The 
abundances of species arising from the TDU are lower, e.g., the 
sum of C-l-N-l-0 (FigurefTOt and the ^--process elements (FigurefTTt. 
The sum of C-l-N-l-0 is determined by the amount of C produced by 
He burning during the TPs, dredged-up to the stellar surface, and 
partially converted to N by HBB. Thus, the truncation of the AGB 
phase naturally leads to a significantly smaller increase in the sum 
of C-l-N-l-0. The higher value of a mlt also has the effect of decreas- 
ing the number of TDU episodes and, in addition, it has the impor- 
tant consequence of altering the thermal structure of the convective 
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Figure 11. Same as Figure [Tol for selected .v-process elements belonging to 
the first (Rb, Sr, Y, and Zr) second (Ba) and third (Pb) i-process peaks. 



envelope. In particular, the temperature at the base of the envelope 
increases. This allows more efficient HBB, which modifies the nu- 
cleosynthetic yields of all the species affected by proton captures, 
for example, the surface O abundance reduces faster (Figure [TOt. 

The i-process elements in our IM-AGB star models are pro- 
duced in the convective regions associated with the TPs. The neu- 
tron source during the TPs is the ^-Ne(Q-,n)-^Mg reaction. The -^Ne 
is synthesised in the TPs by conversion of all the '''N, produced 
by the H-buming shell during the interpulse period, via two a cap- 
tures. An important fraction of the ^^Ne abundance is of primary 
origin since it is formed from '*N which is produced via proton 
captures on the '^C left behind by He burning and carried to the 
envelope by the TDU. The other possible neutron source in AGB 
stars is the '^C(Q',n)"'0 reaction. Formation of the required '^^C is 
usually obtained by allowing some extra mixing of protons at the 
deepest extent of th e convective envelope during each TDU episode 
l lBusso et alJI 19991 ). However, we did not include the "C pocket 
in our models because in IM-AGB stars the formation of the '^C 
pocket appears to be inhibited by p captures occ urring at the hot 
base of the convective envelope during the TDU JGorielv & SiessI 
l2004r) . The surface .^-process abundances increase with each TDU 
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episode as they are successively mixed into the convective envelope 
(see Figure [TTt. The activation of the ^^Ne neutron source depends 
on the temperature at the base of the TP convection zone, which 
increases with each TP. Because the .s-process nuclei accumulate 
in the He intershell and are irradiated during each TP, they expe- 
rience a total neutron flux that increases with the TP number. As 
a result, the j-process surface abundances show an increase ini- 
tially of the elements at the first .?-process peak (Rb, Sr, Y, and Zr), 
then of the elements at the second peak (e.g., Ba), and finally of 
those of the third peak at Pb. The j'-process elements belonging 
to the same neutron magic peak (e.g., Sr, Y and Zr, N = 50) are 
produced together along the .?-process path, according to the lo- 
cal neutron capture rates. In Figure [TT] we show the production of 
heavy elements for two 6 M© AGB models (standard case and aMLT 
enanched-l-Bloecker mass-loss law). For the first model, [Rb/Fe] in- 
creases up to ~1.4 dex, Sr-Y-Zrup to ~ 1.1 dex [Ba/Fe] up to 0.7 
dex and Pb only marginally changes, growing less than 0.2 dex. 
In the second model, the .s-process production is much weaker, 
since all elements under consideration have abundances lower than 
roughly 0.5 dex. To explain variations in only one of these elements 
(eg., Y as reported by VGll) would require very unusual condi- 
tions. The production of Rb with respect to solar is higher than the 
production of Y, due to the high neutron densities resulting f rom 
the activation of the ^-Ne neutron source Jvan Raai et al.ll2012h . 



4.2.3 Comparison with obsen'ations of M4 

In Figures [T2] to |2T| we compare the elemental yields in the form 
of [X/Fe] predicted by the present models with the key abundances 
observed in M4. The "yields" are defined as the integrated elemen- 
tal abundance lost in the AGB winds during the whole life of the 
star. This is the material that presumably contributes to the gas from 
which the SG stars formed. 

- Na and O: In Figures [T2l and [T3l we present the comparison 
between models and observations in the Na-O plane for the 5 Mq 
and 6 Mq models, respectively. It can be seen that the MLT5 model 
covers the range of Na increase and O decrease, while the other 5 
Mq models fail to deplete O sufficiently. On the other hand, most 
of the 6 Mq models deplete enough O, but fail to produce signif- 
icant amounts of Na. The MBR model, computed by altering the 
p-capture reaction rate s that affect the N a abundance (within the 
uncertainties given by llliadis et al.ll2010f) . results in a Na yield ~ 
70% higher than using the recommended rates. This is not enough 
to match the observed variation between the FG and S G stars, how- 
ever, i t should be considered as a conservative choice, llliadis et alj 
(|2010|) provided new estimates for nuclear reaction rates based on 
a Monte Carlo model that allowed the evaluation of the probabil- 
ity distribution function for each rate. They defined as "low" and 
"high" those rates corresponding to a coverage probability of Icr 
- clearly the uncertainties would be larger if one adopted a higher 
coverage probability. Also, it should be noted that while the Monte 
Carlo method takes into account all cuiTent experimental informa- 
tion, it cannot be excluded that different information may arise 
from future experiments, particularly for the ^^Ne(p,7)-'Na and 
^^Na(p,Q')-''Mg rates, which are affected by low-energy resonances 
that are difficult to reach experimentally. 

- Mg and Al: In Figures [14] and [TS] we present the comparison 
between models and observations in the Mg-Al plane for the 5 M© 
and 6 Mg models, respectively. The STD and MLT models produce 
Mg and/or Al abundances too high with respect to the observations, 
either due to dredge-up of ^'-^^Mg produced by a captures on ^^Ne, 



in the case of the 5 Mq models, or due to efficient HBB in case of 
the 6 Mq models. On the other hand, all the BLK and M+B mod- 
els are well within the observations (we remind the reader that the 
extent of the Al spread is debated, as described in detail in Sec- 
tion|4TJ. 

- Na and Y: In Figures [T6] and [T7] present the comparison be- 
tween models and observations in the Na-Y plane for the 5 Mg and 
6 Mg models. The STD and MLT models, which experience a large 
number of TDU episodes, produce Y variations an order of magni- 
tude higher than the spread observed in M4. Therefore, within the 
IM-AGB pollution scenario these models would cause a larger s- 
process enrichment in SG stars compared to FG stars, which is not 
observed according also to our present results. On the other hand, 
the BLK and M+B models, which experience fewer TDU episodes, 
produce Y variations of < 0.3 and 0.2 dex for the 5 Mg and 6 Mg 
models respectively. These models may reproduce better the obser- 
vations, i.e., the missing extra-enrichment in i'-process elements of 
SG stars compared to FG stars in M4 within the uncertainties. A 
similar conclusion can be derived using the Na-Ba plane (Figures 
[l8]and[T9ll. 

- Li: Figures l20l and 127] show the comparison between models 
and observations in the T^^^-Li plane for the 5 Mg and 6 Mq mod- 
els. All our models present a brief phase of Li enrichment at the 
stellar surface, during the first few TPs, with a peak of A(Li)~4 
dex. The final Li yield depends on when the star loses most of its 
mass JTravaglio et aLll200ll ; IVentura & D'Antonall2010l) . The STD 
and MLT models, which experience extended lifetimes, result in 
yields close to zero, because they lose most of their mass after Li 
has been destroyed by HBB. Conversely, the BLK and M-l-B mod- 
els, which have shorter lifetimes, eject a substantial amount of mass 
during their Li-rich phases, producing relatively high yields - close 
to the initial value - particularly when considering the 5 Mg model. 

Overall, none of our STD and MLT models can reproduce 
the observed constancy in the .s-process and Li abundances. Fur- 
thermore, they produce Mg and Al outside the observed variations. 
On the other hand, the BLK and M+B models can reproduce most 
of the observational constraints, although the 5 Mq models do not 
deplete enough O and result in slightly too high i-process abun- 
dances, while the 6 Mq models require some adjustment in the re- 
action rates to explain the Na spread, and they deplete too much 
Li. The Li problem might be rectified by a further slight increase 
in mass-loss rate, for instance using a slightly higher value of r] in 
the Bloecker formula, so that more mass is lost during the Li-rich 
phase, and/or by considering in detail the uncertainties in the nu- 
clear reaction rates involved. 

Finally, we report that in the M+B6 model the [Rb/Fe] in- 
creases of ~ +0.2 dex, slightly higher than that predicted for [Y/Fe]. 
However, this is not a problem since such a weak overall enrich- 
ment of Rb and Y in these IM-AGB models as polluters would 
likely not be observable in SG stars. For the same model, the 
change in the sum of C+N+O is — hO.3 dex, and the He enrich- 
ment is 6Y = +0.09. All the models produce [N/Fe] ^1.5-2 (the 
observed variation is up to =:I dex; e.g., VGl 1), but none of them 
deplete C to the observed level of [C/Fe]^-0.5 (e.g., VGll). The 
lowest [C/Fe] obtained is =: using the M+B6 model. 



4.2.4 Comparison with other observational constraints 

To ascertain if our IM-AGB models computed including efficient 
HBB and strong mass loss are a realistic solution for the anoma- 
lies observed in M4 we can, in principle, check if these models are 
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Figure 13. Same as Figure [Tll except for the 6 M© models. 
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have offset the initial Al composition used in the models to the value 
observed in the FG stars of M4. 
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Figure 15. Same as Figure fMl but for the 6 M© models. 



consistent with other observational constraints. A problem how- 
ever is that there are very few direct observations of abundances 
in IM-AGB stars, and these are for stars in the Galaxy, LMC, or 
SMC, i.e., of meta llicities higher than M4. One study that we are 
aware of is that of iMcSavenev et alj OOOTr) . who obtained spectra 
of a sample of LMC and SMC AGB stars. The se vere difficulties in 
perfor ming these observations are clear from the iMcSavenev et alj 
( 120071) study, where they attempted to measure abundances in the 
cool AGB atmospher es (which are also dynamic/pulsating, see e.g., 
iLebzelter et al.l2010r) of 7 stars but were only able to obtain results 
for two. The results from their two stars support an increase in omlt 
at least as high as we have used in our MLT models and also the 
phenomena of HBB and TD U as seen in the m odels. Other studies 
of IM-AGB stars are those of IPlez et al"* ('1993') for t he Smal l Mag- 
ellanic Cloud (SMC), and those of iGarcia-Hemandez et al.l (2006, 



2007, 2009) discussed below, for the Galaxy, the SMC and the 
Large Magellanic Cloud (LMC). More accessible are planetary 
nebulae (PNe) of Type L The material surrounding the compact 
central stars in these objects is believed to consist of the mate- 
rial ejec ted by IM-AGB stars. They are also located in the G alac- 
tic disk JStangheUini e t al."2006 l:ISterling & Dinersteinl2008h . or in 
the LMC and SMC fcaler & Jacob'^ll99Cl) . Again, all these ob- 
jects have metallicities higher than M4, which may affect the un- 
certain parameters related to the evolution of these stars. Keeping 
this caveat in mind, we can examine the available constraints. 

Planetary nebulae (PNe) of Type I are believed to be IM stars 
that have completed the AGB phase, have little or no stellar en- 
velope left, and are evolving blue-wards at constant luminosity in 
the CMD. PNe of Type I show the clear signature of HBB in their 
He/H and N/0 ratios (e.g.. lKaler & Jacobvlll990l : IStanghellini et alj 
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Figure 17. Same as Figure [T6l but for the 6 Mq models 
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Figure 19. Same as Figure [Tsl but for the 6 M© models. 



120061) and no enhancements of the .s-process elements Se and Kr, 
which are close to the first i'-process peak. This is well explained 
by m odels of IM-AGB stars of solar metallicity JKarakas et alj 
1 20091) . However, this simple picture has been somewhat compli- 
cated by the detection of extremely high abundances o f Rb in 
OH/IR stars in t he Galaxy JGarcia-H emandez et al.ll2006h . and in 
the SMC/LMC JGarcia-Hemandez et al...200!^) . which represents 
the first proof of the activation of the -^Ne neutron source in IM- 
AGB stars. The IM-AGB nature of these stars was determined on 
the basis of the velocities of their OH masers, their location in 
the Galactic plane (as indicative of belonging to a young popula- 
tion), and the presence of Li produced by HBB, at least in some 
of the objects. The very high Rb overabundances, with [Rb/Fe] 
up to -1-2.5 dex for Galactic stars and up to 4-5 dex for stars in 
the Magellanic Clouds are however not matched by standard IM- 



AGB mod els tearcia-Hemandez et al."2009Vvan Raai et al."2012|). 
IKarakas et al. ( 2012) showed that if the stellar lifetime is extended 
by lowering the mass-loss rate, so that a larger number of TDUs is 
allowed, and the rate of t he "^Ne(Q',n)^^Mg r eaction is taken from 
the NACRE compilation jAngulo et alJI 19991 ). then it is possible to 
reach [Rb/Fe] ^ 1.4, close to the average [Rb/Fe] observed in the 
Galactic OH/IR stars. Clearly, we have here a contradictory situa- 
tion since to match these observations the .y-process nucleosynthe- 
sis in IM-AGB stars needs to be enhanced, while to match the con- 
stancy of the .s-process elements in M4, including Rb, the .s-process 
nucleosynthesis in IM-AGB stars needs to be suppressed. 

We note however that there are still several problems related 
to the interpretation of the high Rb abundances in IM-AGB: 



(i) as discus s ed in detail bv lGarcia-Hemandez et alj ( 120091) and 
Ivan Raai et alj ( l2012n . serious problems are present in current 
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Figure 20. Same as Figure [T2l but for the T^ff-Li plane. The thick, soUd 
line is the initial abundance of Li, the other hues represent the predicted 
yields as indicated by the labels. Observational data are the dwarf stars 
discussed in Section l4~l1 
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Figure 21. Same as Fig |20l for the 6 M© models. 



model atmospheres of luminous AGB stars since these models 
are performed in ID and do not include important effects such as 
the presence of a circumstellar dust envelope and dust formation, 
which may lead to systematic uncertainties. More realistic model 
atmospheres for intermediate-mass 0-rich AGB stars (e.g., the in- 
clusion of a circums tellar dust envelope a nd 3D hydrodynamical 
simulations, see e.g.. lLebzelter et al.ll2010t) as well as NLTE calcu- 
lations need to be developed; 

(ii) the IM-AGB models that predict high Rb enhancements also 
necessarily predict some Zr enh ancements, which are not obser ved. 
These stars show [Zr/Fe] < 0.5 JGarcia-Hemandez et al.ll2007l) : 

(iii) the NACRE -^Ne(Q',n)-^Mg rate i s probably too high when 
considering the more recent experiments (IJaeger et al.bOOlh : 

(iv) Li abundances in the observed stars present a large range, 
from depletion to enhancements, however (as also discussed above 
in relation to the present models), Li enhancements can only be 
seen early during the AGB evolution, while high Rb enhancements 
are obtained only after many TPs. 

Another possibility is that the TDU efficiency and the mass- 
loss rate are very sensitive to the initial stellar mass and the metal- 
licity. For example, more m assive IM-AGBs than those considered 
here and Su per-AGB stars dPumo et al.ll2008l iDohertv et al.ll201^ 
ISiesj 12010!) might experience less efficient TDU and/or higher 
mass-loss rates and be responsible for the pollution of GCs. Obvi- 
ously, more models of IM-AGB stars and observations of both Rb- 
rich stars and Type I PNe will help us to understand these issues. 
For the time being it seems difficult to use the high-Rb observations 
as a strong constraint to extrapolate information to IM-AGB stars 
in GCs and viceversa. 

A further difficulty when comparing to other observational 
constraints is the possibility that IM-AGB stars could evolve dif- 
ferently in clusters and in the field. It is known that B stars (which 
eventually evolve to be IM-AGB stars) in open clusters typically 
rotate more rapidly than B stars in the field (though this may be 
due to their being in a less evolved phase of their evolution, e.g., 
iHuang et alJl201CT) . Fast rotation would have a strong effect on the 
stellar evolution prior to the AGB phase possibly producing too 
high C+N -hO yields to be comp atible with the observations, as 
shown by iDecressin et alj ( 120091) , though the effect of magnetic 



fields was not included in these models. Mor eover, binary inter- 
action is important in shaping stellar yields ( llzzard et al.l l200q , 
IVanbeveren et al.ll2012h and its effects would be different if the bi- 
nary distribution was different for field and cluster stars. 



5 SUMMARY AND CONCLUDING REMARKS 

In summary we have considered four observational constraints and 
derived the following results for M4: 

(i) In the O-Na plane, the MLT5 model provides a good fit to the 
observed spread. Also the STD6 and the M+B6 models may result 
in a good fit, provided that the rates of the nuclear reaction that pro- 
duce and destroy Na during HB B are varied somewh at beyond the 
'high' and 'low' values given bv llliadis et alj ( 120101) . These varia- 
tions are within current possibilities and need to be investigated. 

(ii) In the Mg-Al plane, all the STD and MLT models produce 
variations in Mg and/or Al outside the observed spread, while all 
the BLK and M+B models are within the observations. 

(iii) As in item (ii), in the Na-Y and Na-Ba planes the STD and 
MLT models produce variations outside the observed spread, while 
the BLK and M+B models can explain the constancy of .s-process 
abundances between generations, particularly when considering the 

6 Mq models. 

(iv) As in items (ii) and (iii), when considering Li the STD and 
MLT models result in no Li production, while the BLK and M-l-B 
models can explain the constancy of the Li abundance between gen- 
erations, particularly when considering the 5 Mq models. 

(v) All the models result in production of He beyond the ob- 
served variations, which indicates the need of some dilution with 
primitive material. The change in the sum of C-l-N-l-O may be com- 
parable to the observations, within the uncertainties, only if we con- 
sider the BLK and M+B models. The observed depletion of C by a 
factor ~ 2 is not obtained by any of the models. 

Although none of these models can simultaneously match all 
the constraints listed above, we believe that our exploration of the 
model uncertainties shows that IM-AGB stars cannot be ruled out 
as potential polluters for the M4 SG population since reasonable 
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variations in input physics can still provide a reasonable agree- 
ment between the theoretical yields and observations. More work 
is clearly needed to improve the input physics of IM-AGB stars. 
We note, however, that the fact that SG stars probably formed from 
the ejecta of a range of masses (not a single mass) and that there 
must be a certain dilution b etween pristine and polluted material 
(see e.g., lGratton et alj2012l) . do not allow us to simultaneously re- 
produce all the observational constraints any better with the present 
models. 

We have c o nfirmed the previous finding by 
IVentura & D'Antonal ( l2Q05cl : 2010) that IM-AGB models 
computed including efficient HBB and strong mass loss (e.g., the 
Bloecker mass-loss rate) can provide a match to the observations 
of the light elements. Furthermore, we have shown that these same 
models can also be consistent with the constraint that there are 
no variations in the .?-process elements in M4 between FG and 
SG. When comparing to observational constraints other than GCs, 
however, it appears difficult to establish a self-consistent scenario 
that could explain all the current observations related to IM-AGB 
stars. Also the question remains of why IM-AGB stars in GCs 
should experience a stronger mass loss than what is observe d in 
the G alaxy and in the Magellanic Clouds ( Vassiliadis & Wood 
Il993h . In any case, pursuing such comparison offers us the future 
opportunity of better understanding all of the current issues related 
to nucleosynthesis in IM-AGB stars via more refined models and 
observations. 

In particular, AGB i-process models similar to those we have 
explored here are needed for a more extended range of masses, i.e., 
between 3 Mq and 9 Mq, and for different choices of the mass- 
loss ra te and efficiency of convection. The models of iLugaro et alj 
( 120 12n have shown that as the stellar mass decreases the effect 
of the TDU (producing C, F, and i-process elements) gradually 
becomes predominant over the effect of HBB (depleting O, and 
F, and producing Li). AGB models of slightly lower mass than 
those presented here have been invoked as a possi ble explanation 
for variations in the .5-p rocess elements in M22 JRoederer et al.l 
I2OI iLlMarino et al.ll2009l). as well as for their correlation with flu- 
orine JD'Orazi et alj|2013n . A large grid of AGB mo dels and real- 
istic model s of the gas dynamics in the forming GC ( IKrause et alj 
l20I2ll20I3h will be fundamental to ascertain at which mass we can 
see the raising of the .s-process and fluorine abundances, and to 
derive information on the timescales of formation of the different 
populations in GCs. 
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